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Abstract. Application of binary linear programming can be found in many fields such as scheduling, production 
planning, networking and etc. Unlike linear programming problems, binary integer linear programming problem 
is more difficult to solve since its variables must be a binary number zero or one. Thus, solving this problem 
efficiently becomes an interesting issue. This paper proposes a combination of Balas additive method and 
Gomory cutting plane method. Balas additive method is a well-known method in binary linear programming. It 
systematically enumerates subset of the possible binary solution to ensure all possible solutions are examined. 
In this paper, a drawback of this method which is requiring long iterations is fixed. It uses the Gomory cutting 
plane method to eliminate non-integer partial solution. The Gomory cutting plane method generates an 
additional constraint to eliminate a non-integer solution. In general, this additional constraint increases the 
problem complexity. Therefore, combining Balas algorithm and Gomory cutting plane method can solve binary 
linear problem faster without facing more complex mathematical model. Substituting Gomory cutting plane to 
the Balas algorithm is conducted by performing the cut based on condition in Gomory algorithm towards the 
enumeration process conducted by Balas algorithm. The computational result using a test function proves that 
the proposed substituting strategy can solve binary linear programming efficiently which is shown by the 


number of iterations required by this method. 


1. Introduction 


Linear programming (LP) has been applied in many different fields. One of application of 
linear programming can be founded in a problem requiring a “yes-no’ decision. Herein, a 
binary zero-one LP can be applied to solve this kind of problem. Due to its importance, many 
algorithms have been proposed to solve a zero-one LP. One of them is Balas algorithm 
proposed by Balas [1]. This method has been proven can solve zero-one LP well, yet it 
requires a long iterations[2]. On the other hand, A Gomory cutting plane method was 


proposed for solving an integer LP[3]. This method does not require too long iterations as 


2010 Mathematics Subject Classification: 90C05. 
Key words and phrases: Binary linear programming, Balas algorithm, Gomory cutting plane. 
©2014 Science Asia 
t711 


2/11 
ROFLIN, OCTARINA, BANGUN, ATMANEGARA AND ZULVIA 


proceed by Balas algorithm. The approach taken in the Cutting Plane technique is to create an 
additional constraint that cuts feasible space of the linear programming relaxation to 
eliminate a solution that is not integer. But this method has the disadvantage such a solution 
obtained was shaped shards or fractional integer which does not require complex simplex 


iteration[4]. Unfortunately, this method does not obtain an absolute zero-one solution. 


This paper aims to overcome the drawback of Balas algorithm by embedding it with Gomory 
cutting plane method. The main concept of the proposed method is eliminating the 
backtracking procedure in Balas algorithm and replacing it with Gomory cutting plane. 

The rest of this paper is organized as follows. Section 2 briefly reviews the literature study of 
Balas algorithm and Gomory cutting plane method. In Section 3, the proposed combination of 
Balas and Gomory algorithm is discussed. Furthermore, two numerical examples are given in 


Section 4. Finally, a concluding remark is given in Section 5. 


2. Literature Review 

This section discusses some basic theories of linear programming, Balas algorithm and 
Gomory cutting plane which are applied in this paper. 

2.1 Binary linear programming 

Linear programming is a tool for solving optimization problem. It aims to minimize or 
maximize the objective function with respect to some constraints. In the real-world problem, 
optimization problems usually involves a “yes-no” decision. These problems are modeled as 
binary linear programming or zero-one linear programming. The general form of binary 
linear programming model is as follows. 

Objective function: 


Maximize or Minimize 


f(x)=C'X (1) 
Subject to: 

AX+Y=B (2) 
ION (3) 


Y>0 (4) 
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foul ea feat a a, ain | 
x c b a a a 
wherex =}, C=) })20,B=);),and A = R a 
| | | | el | | 
[xa] |c, | [ip] | ni d n2 Ginn | 


2.2 Balas Algorithm 

Balas algorithm was proposed by Egon Balas [1]. This method works by successively 

assigning the value 1 to a certain variables. Notations used in Balas algorithm are defined as 

follows 

Jy = Set of decision variables at iteration & Itisa partial solution at iteration X. 

R = Set of decision variables in objective function 

Jk = Set of decision variables which are not included in J,. All decision variables are 
assumed the be zero unless it has been defined as one in advance.J; = R — Jp. 

P = Set of decision variables which can repair the final solution. 

IX = Number of violated constraints when variable x, = 1 at iteration &. 

fx = Objective value at iteration &. 


f° = Final objective value 


yj = Slack variable at iteration k. 


The Balas algorithm is as follows. 

STEP 1. Initialization. 

Set k= 0, J) = Ø, YY= B with at less one component from B< 0, f* =o, fa = 0, R= 
{Xi X>, Xn} and Jo =R. 

STEP 2. Determine variable in J, which can be moved to Jẹ. 

Choose one set of decision variables in Pwhich has objective value better than Jx Eliminate 


variables which have negative value. 


yy ea ee (5) 


Variable yj is feasible if ajs < 0. The x, € J. With aj; = 0 can be eliminated and inserted to 


the set Nx Eliminate all variables causing the current objective value f, higher than f*. Let 


the objective function: 


ae corse (6) 


jeJ, 
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If x, E J,, is set as 1, the objective function becomes (fy + c,). Thus x,which makes the 
objective value greater than f*cannot be included in solution. Let M, = {x,|x, E L, L = fk + 
Cs => f*}, then: 


P.=J,-(N,U¥M,)=R-J,-(N,UM,) (7) 


If P, = Ø, then go to STEP 4. Otherwise, go to the next step. 


Given constraints: 
A,X, +a,,X,+...+ a,x, + Y, = Ose NG <0 (8) 
which can be written as: 
e ay, SO (9) 
x EJ; 
By considering only variables in P,, we get: 


D ap SY; VY; <0 (10) 
x EP, 


Change variable,x, € Px, to 1. The variable x, is chosen based on Eqs. (11). 


x, = argmin (r; = max(7‘)} (11) 
where, 
i vs ans aa ee P: (12) 


j=l 


Defines new solution J;44. 


STEP 3. Update y;**?: 


y 


ee ee (13) 


J 


Update objective function: 


Pa re (14) 
If all aig > 0 ,then f* = fk+1, and go to STEP 4. Otherwise, back to STEP 2. 
STEP 4. Backtracking and determining optimal solution. 


The backtracking process is conducted towards partial solution Ję. It is started by assuming 


all variables in J, are one. Then, successively assigning zero to certain variables. Herein, 2” 
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possibilities are evaluated where n is number of variables in Jẹ. Then, a combination which 


gives the smallest value is the final objective value. 


2.3 Gomory Cutting Plane 

The Gomory cutting plane is an algorithm for obtaining integer solutions of the linear 
programming [3]. This algorithm works by examine a solution of the linear programming 
obtained by simplex method. If the solution is not in integers, a new constraint which can cut 
the search space so that non integer solution can be eliminated. The procedure of generating 
new constraints is as the following. 

STEP 1. Given an optimal simplex tableau. 


If x; < 0 , then 


Z alw, =B, (15) 


when|a/ | < al andw; = 0, equation (15) can be transformed to: 


5 alw, Sf. (16) 
and 
5 alw, +S =P. (17) 


where S; is a slack variable. Assumed f; is non-integer, then|a? | + fij = al and |i] + fi = 


Piwhere 0 < fi; < land 0 < fi < 1. The additional constraint is the differences between Eqs 
(15) and (17) which can be formulized as follows. 


ree a (18) 


Jal 


3. Methodology 

The idea of substituting Gomory cutting plane to the Balas algorithm is proposed to accelerate 
the procedure of Balas algorithm. As discussed in Section 2, Balas algorithm involves a 
backtracking procedure which takes a long iteration. Thus, we perform Gomory cutting plane 
method to replace the backtracking procedure. Herein, Gomory cutting plane condition is 
applied to determine which variable should be set as 1. Detail steps of the proposed Balas and 


Gomory combination is given as follows. 


6/11 
ROFLIN, OCTARINA, BANGUN, ATMANEGARA AND ZULVIA 


STEP 1. Initialization 

Set k = 0. Transform the LP to a standard form. 

STEP 2. Solve the standard form LP using dual simplex method to obtain the optimal dual 
simplex tableau. 

STEP 3. If all decision variables in the optimal solution are integer, then stop. Otherwise, 
generate a new constraint based on Gomory procedure. The new constraint is derived from 
the basis variable which has the biggest non-integer solution. Solve the new LP with dual 
simplex method and obtain new optimal tableau. If the new solution is integer, then stop and 
go to STEP 4. Otherwise, back to STEP 3. 

STEP 4. Substitute solution the original LP form to get the final optimal solution. 


4. Numerical Analysis 

In order to analyze the performance of the proposed method, a numerical analysis involving a 
binary LPs is conducted. Herein, a binary LP is solved using Balas algorithm, Gomory cutting 
plane method, and the proposed method. 

Problem 


Objective function: 


MIN f =-5x,+7x, +10x, —3x, +x, (19) 
Subject to: 

x, +3x,—-5x,+%x,+4-, < 0 (20) 
2x +6x, -3x +2x,+2x, 24 (21) 
E ERES E (22) 
veloi Vis iss (23) 


4.1 Procedure the proposed substituting Balas and Gomory 
Solve the LP using the proposed Balas and Gomory algorithm as follows: 
Iteration 1. 


STEP 1. Transform LP to standard form and obtained new form as follows: 
MIN f =—5x,+7x, +10x, —3x, + x, (24) 


Subject to: 
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—~x,+3x,—-5x,-x,+4x,+y,=-2 (25) 
2x, —-6x,+3x,+2x,-2x,+y, =0 (26) 
x,—-2x,+x,+x,+y,=-1 (27) 
x, e {0,1 Vi=1,...,5 (28) 
Y Y, Y, 20 (29) 


STEP 2. Solve the standard form using dual simplex method and obtain the optimal tableau as 


follows. 


| Baic | ox, |n| n] u | x | y | y: | ya | Solution | 


| 
“Eke 
| 
‘D| 
| | 
MO! wkol m 
| 
ail H/o] bo 
H 


= 
ONO0EIE 
9 
7 
Pa Pepa Pa 


Figure 1. Optimal tableau 1 


1 es 
le 





Since the current solution contains non-integer variables, a new constraint must be generated. 
According to the Gomory cutting plane rule, the x, is used for the new constraint. The new 
constraint generating procedure is as follows: 

Change negative coefficient of all variables in row x3 into positive value. 

x, > -2/9=7/9-1 

x, > -8/9 =1/9-1 

x, > -4/9 =5/9-1 (30) 

y, > -1/9 =8/9-1 

y, > -2/3 =1/3-1 


Define new coefficient for all variables: x, =7/9;x, =1/9;x, = 5/9; y, =8/9; y, =1/3 


Substitute to equation x3. 


7 1 5 8 1 2 
M a a a a (31) 
9 9 9 9 3 3 


Define the constraint as follows: 
7 1 5 8 1 2 
So Gg a et a a (32) 


Add the new constraint to the optimal tableau and solve using simplex method. The result is 


shown as follows. 
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Bp o eee ee fon ee 





MOH ecee 
4 


Figure 2. Optimal tableau 2 


From the current optimal tableau, the optimal solution still contains non-integer variable. 
Thus, return to STEP 3. 

Iteration 2. 

STEP 3. Generate a new constraint. 

From current optimal tableau, xand y, are variables which have the biggest integer 
solution. However, since coefficient s4 in row y, is less than -1, the new constraint is 
generated fromx3. By following the same procedure in iteration 1, we get the new constraint 


as follow. 


33 
a (33) 
+) ele |S 
os) 
— 
“a 
ET eee ee ee (34) 

8 8 8 8 8 4 


Add the new constraint to the current LP and solve it using Simplex method. The new optimal 


tableau is as follows. 


x| xa| x, | xs |y, |y | y [| S| S, | Solution | 


BIBOORHOOMOEs 
EInCoesHOoRoeee 


Figure 3. Optimal tableau 3 





Since the optimal solution still involves non-integer solution, go back to STEP 3. 


Iteration 3. 
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STEP 4. Generate new constraint. 
The optimal tableau shows the highest integer values are given by xand s,. However, s, in 


row sı has value less than -1. Thus, x3 is chosen for the new constraint as follows. 


x ay pe esla + PAN 

a a 

e[t- Ss, (35) 
(. 6) 

tae 

$,-—2,-24,-—9,- 28, =-— (36) 


Add the new constraint to the LP and solve it using Simplex method. The optimal tableau is as 


follows 


se [afn a| u te tell me tt | Se | Soton 





S44 
= 
5 


Figure 4. Optimal tableau 4 


The optimal tableau shows that the solution is non-integer. However, there is one variable has 
negative value. Thus, recalculating coefficient f and y, is performed instead of generating a 


new constraint. The ratio between coefficient f and y, is given as follows. 
PPE EE ES 


2 a 2 4 
Yı -1 —3- 1 —3- —1-| 1- 
5 5 5 5 
ss 5 5 
Ratio | 5 4— 2 1— 
6 7 


Figure 5. Ratio between coefficients f and yı 





Choose variable with the smallest ratio to replace y, is basis row. Recalculate the optimal 


tableau using Simplex method. The result is as follows. 
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2 4 5 5 2 
| f (zleje faol o|a3l0| -7| 0 e=3] "| 
waa aea 
Xa -- 1 | -1- = -- 
pe | 5 [tfel ag fe] -pfel ago feist a 


| | 5 fefe| 25 [ols lel ste ts fas] 
» [szele] 75 [elez sze fol] > | 
3 6 3 4 a, 
s | = 4- = 4- Ao 2 
ps | 5 Jefol-25]s] 5 fol-s}e jel se] o 
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Figure 6. Final optimal tableau 


This optimal tableau shows that the current solution is integer solution. Thus stop the 
iteration and go to the next step. 
STEP 4. Substitute the optimal solution to the original LP form to get the final optimal 


solution. 


x, =x, =S8,=l,y, =3,8, =2,x, =0 (37) 


Soe Y 
x =x, =x, =x, =l (38) 
x. =0 


4.2 Discussion 

The same problem is also solved using Balas and Gomory algorithm. Since the aims of 
combining Balas and Gomory algorithm is to reduce the computational time in terms number 
of iterations, the evaluation is conducted by comparing number of iterations required by 
Balas, Gomory and the proposed method. The result is mentioned in Table 1. This result 
proves that the proposed Balas and Gomory algorithm can solve the binary LP better than 
Balas and Gomory algorithm. It is not only can obtain better result, but also only requires less 
iteration. 


Table 1. Comparison between Balas, Gomory and the proposed method 


Method Number of Objective 
iterations 


Balas Algorithm 5 iterations: 32 


steps 
1 iteration: 2 


steps 


The proposed method 
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5. Conclusion 

Balas algorithm is a method for solving binary LP. Although it has been proven can solve the 
binary LP well, this method has a drawback which is requiring a long iterations. On the other 
hand, Gomory method is a method for solving integer LP. Unlike Balas algorithm, this method 
only takes a few iterations to find an integer solution. 

This paper aims to propose a method for solving binary LP with less iteration. The proposed 
method combines Balas algorithm with Gomory cutting plane method. In this method, the 
backtracking procedure in Balas algorithm is replaced with Gomory cutting plane method. 
Backtracking procedure is a part of Balas algorithm which requires 2” evaluation where n 
is number of variables. The proposed method avoids this procedure since it requires high 
computational. Herein, Gomory condition is applied to generate an additional constraint 
which reduces the search space so that it can find a binary integer solution faster. 
Furthermore, the proposed method is validated using a binary LP. Numerical computational 
result proves that the proposed substitution Balas and Gomory algorithm can solve binary LP 


better and faster than Balas and Gomory algorithm. 
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